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A key ingredient in social contagion dynamics is reinforcement, as adopting a certain social behavior requires 
verification of its credibility and legitimacy. Memory of non-redundant information plays an important role in 
reinforcement, which so far has eluded theoretical analysis. We first propose a general social contagion model 
with reinforcement derived from non-redundant information memory. Then, we develop a unified edge-based 
compartmental theory to analyze this model, and a remarkable agreement with numerics is obtained on some 
specific models. Using a spreading threshold model as a specific example to understand the memory effect, 
in which each individual adopts a social behavior only when the cumulative pieces of information that the 
individual received from his/her neighbors exceeds an adoption threshold. Through analysis and numerical 
simulations, we find that the memory characteristic markedly affects the dynamics as quantified by the final 
adoption size. Strikingly, we uncover a transition phenomenon in which the dependence of the final adoption 
size on some key parameters, such as the transmission probability, can change from being discontinuous to 
being continuous. The transition can be triggered by proper parameters and structural perturbations to the 
system, such as decreasing individuals’ adoption threshold, increasing initial seed size, or enhancing the network 
heterogeneity. 

PACS numbers: 89.75.Hc, 87.19.X-, 87.23.Ge 


I. INTRODUCTION 

Due to technological advances social networks are playing 
an ever increasing role in the modern society. In a social net¬ 
work, nodes are individuals of the population while links rep¬ 
resent the social ties or relations among individuals Et In 
recent years, there is a growing interest in investigating the 
phenomenon of behavior spreading on social networks, where 
the behaviors range from adoption of an innovation El] and 
healthy activities [ 3j] to micro finance 0], This is essentially 
the problem of social contagion. Ample experimental and 
theoretical results indicated that, unlike biological contagions 
in which successive contacts result in contagion with inde¬ 
pendent probabilities, in a social contagion the probability of 
infection depends on previous contacts. This is equivalent to 
social affirmation or reinforcement effect, since multiple con¬ 
firmation of the credibility and legitimacy of the behavior is 
always sought ffl. For an individual, who had two friends 
adopting a particular behavior before a given time and whose 
third friend newly adopts the behavior, whether he/she adopts 
this behavior will take the three friends into account. 

An early mathematical model to describe the dynamics of 
social contagions is the threshold model IIIOL lllll based on 
Markovian process without memory, in which adoption of be¬ 
haviors depends only on the states of the current active neigh¬ 
bors (i.e., individuals who have adopted the behavior), and an 
individual adopts a behavior only when the current number or 
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the fraction of his/her active neighbors is equal to or exceeds 
some adoption threshold. Analytically, the fraction of individ¬ 
uals adopting the behavior eventually, can be predicted using 
the percolation theory ifTTIl for situations where the initial seed 
size is vanishingly small. One result is that, for a fixed thresh¬ 
old, as the mean degree is increased, the final size tends to 
grow continuously and then decrease discontinuously. As the 
degree distribution becomes more heterogeneous, the network 
is less vulnerable to social contagions, in sharp contrast to the 
dynamics of epidemic spreading EMI. Previous research 
also revealed that, within the threshold model, factors such 
as the initial seed size M, clustering coefficient Ezth com¬ 
munity structure EMm multiplexity ||20l-l22ll. and temporal 
networks i23il24l l all can affect the social contagion process. 

In real situations of social contagions, memory typically 
plays an important role in the adoption and reinforcement of 
behaviors, which includes full El] or partial [H memory of 
the cumulative behavioral information (behavioral informa¬ 
tion can be referred as information for short) that individuals 
received from their neighbors. This memory effect makes the 
dynamics of social contagions have non-Markovian charac¬ 
teristic. To account for the memory effect, sophisticated non- 
Markovian models were proposed IBLHi, 6i, 171 125142711 . In some 
models, it was predicted that the final adoption size will grow 
discontinuously 7j, ill , if the adoption probability for an 
individual who receives more than one piece of information is 
two times larger than the adoption probability for individuals 
getting only one piece of information. In general, the memory 
of cumulative information about the particular social behavior 
can come from redundant SH or non-redundant El infor¬ 
mation transmission, where the former allows a pair of indi¬ 
viduals to transmit information successfully more than once 
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but for the latter, repetitive transmission is forbidden. In some 
social contagion processes such as risk migration and use of 
unproven technologies @], transmitting redundant informa¬ 
tion between the same pair of individuals is unnecessary, since 
each neighbor can guarantee the credibility and legitimacy of 
the behavior but only to certain extent [3|]. However, such 
non-redundant information transmission characteristic of so¬ 
cial reinforcement have essentially been neglected in previous 
studies IHi- 

A systematic study to understand the effects of non- 
redundant information memory on social contagion dynamics 
is thus called for. A general model needs to include differ¬ 
ent situations of behavior adoption such as the dependence 
of the adoption probability on non-redundant information [ 3i] 
or even on the structure diversity of such information 12811 . 
Due to the non-Markovian nature of the memory character¬ 
istic, to develop a general theory is challenging. Some ap¬ 
proximate approaches were devised such as those based on 
mean field analysis 0], percolation theory [[251], and renewal 
process [295, 13011 ■ Since the non-Markovian property induces 
strong dynamical correlations between any two connected in¬ 
dividuals, analytic predictions from these approaches tend to 
deviate significantly from results from direct numerical sim¬ 
ulations, especially when the underlying network is strongly 
structurally heterogeneous & 

In this paper, we articulate a general social contagion 
model with social reinforcement derived from memory of 
non-redundant information to address the general question of 
how behaviors spread on networks in a more systematic and 
complete way. In order to understand, quantitatively, the ef¬ 
fects of this kind of memory characteristic on social conta¬ 
gion dynamics, we develop a unified edge-based compart- 
mental theory. We base our study on the spreading thresh¬ 
old model, focusing on the final behavior adoption size and 
its dependence on the transmission probability under different 
dynamical and topological parameter settings. We find that 
the memory characteristic generally have a strong effect on 
the final adoption size. Surprisingly, we uncover a crossover 
between discontinuous and continuous variations in the final 
adoption size. More specifically, the crossover phenomenon 
can be induced by decreasing individuals’ adoption threshold, 
increasing the initial seed size or enhancing the structural het¬ 
erogeneity of the network. Our theoretical predictions agree 
well with results from numerical simulations. We further gen¬ 
eralize our theory to treat distinct social contagion models and 
network structures. 

In Sec. HU we describe our general social contagion model 
with reinforcement derived from memory of non-redundant 
information on complex networks. In Sec. [Till we detail our 
edge-based compartmental theory and analysis. In Sec. lIVI we 
present results from extensive numerical computations to val¬ 
idate our theory. In Sec. [VJ we extend our theoretical frame¬ 
work to analyze alternative social contagion models, demon¬ 
strating the generality of our theory. In Sec. [VT] we present 
conclusions and discussions. 


II. A GENERAL SOCIAL CONTAGION MODEL 

Our goal is to construct a general stochastic model for 
social contagion dynamics, taking into account social rein¬ 
forcement through non-redundant information memory char¬ 
acteristic. In this model, information refers to the behav¬ 
ioral information. The non-redundant information memory 
has two features: (1) non-redundant information transmission, 
i.e., repetitive information transmission on every edge is for¬ 
bidden, and thus also can be called as single-transmission ; 
(2) every individual can remember the cumulative pieces of 
non-redundant information that the individual received from 
his/her neighbors, which makes the contagion processes be 
non-Markovian. 

Concretely, we consider a configuration network 
model [32] of size N and degree distribution P(k), 
where nodes in the network represent individuals. There is 
no degree-degree correlations when the network is very large 
and sparse. At any time, each individual can exist in one of 
the three different states: susceptible, adopted, or recovered. 
In the susceptible state, an individual does not adopt the 
social behavior. In the adopted state, an individual adopts the 
behavior and transmits the behavioral information to his/her 
neighbors. In the recovered state, an individual loses interest 
in the behavior and will not spread the information further. 
This is thus a susceptible-adopted-recovered (SAR) model. 
Although this proposed model has similar state definitions 
with the epidemiology susceptible - infected - recovered 
(SIR) model iTH, the non-markovian characteristic is absent 
in the SIR model. 

To initiate a social contagion, a fraction p (J of individuals 
are uniformly randomly chosen to be in the adopted state and 
the remaining majority of the individuals are in the suscep¬ 
tible state. At each time step, behavioral information propa¬ 
gates from each adopted individual to each neighbor indepen¬ 
dently with transmission probability A, a key parameter of the 
underlying dynamical process. We assume an edge that has 
transmitted the information successfully will never transmit 
the same information again, i.e., non-redundant information 
transmission. The non-redundant information spreading pro¬ 
cess is illustrated schematically in Fig. 02 a). Based on this 
setting, we introduce the memory effect of non-redundant in¬ 
formation in social reinforcement. In particular, assume that a 
susceptible individual u of degree k already has rn - I pieces 
of information from distinct neighbors. Once u is success¬ 
fully informed of the social behavior by one of his/her adopted 
neighbors, denoted as v, the cumulative number of pieces of 
information that u has will increase by 1. With the m cumu¬ 
lative pieces of information up to now (i.e., after exposing to 
to pieces of non-redundant information), the probability that 
the individual will be in the adopted state is n(k,m). Note 
that u may subsequently get more than one pieces of infor¬ 
mation successfully in this time step, thus, he/she will try to 
adopt the behavior when he/she gets every new piece of infor¬ 
mation. In this case, if u gets the (n + l)th new information 
in this step, he/she will adopt the behavior with probability 
7r (k,m + n). An illustration of the behavior adoption pro¬ 
cess is presented in Fig.QJb). Since 7 r{k,m) < 1 in general. 
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FIG. 1. (Color online) Illustration of the susceptible-adopted- 
recovered (SAR) model on complex networks, (a) At time t, the 
adopted individual 5 tries to transmit the behavioral information (or 
simply information for short) to each susceptible neighbor individ¬ 
ual independently with probability A. Note that individual 5 can not 
transmit the information to individual 4, since he/she has transmitted 
the information to individual 4 successfully before time t. That is to 
say, the susceptible individual can only get the non-redundant infor¬ 
mation from his/her neighbors. The solid blue lines denote that the 
information has not transmitted through them successfully, and the 
red dished line denotes that the information has transmitted through 
it previously, (b) Assuming that individual 1 has received a new 
piece of information at time t, whether individual 1 adopts the behav¬ 
ior is determined by the m cumulative pieces of information he/she 
ever received from neighbors. The value of m can be expressed as 
m = md, where rrid is the pieces of information that indi¬ 

vidual 1 received at time d. In such a situation, individual 1 has to 
remember the pieces of non-redundant information he/she received 
from neighbors before time t. Thus, the so called non-redundant in¬ 
formation memory is induced. Individual 1 becomes adopted with 
probability 7r(fc,m), where k is the degree of individual 1; other¬ 
wise, individual 1 remains in the susceptible state. 


multiple information transmission is necessary for u to move 
into the adopted state, thereby incorporating the memory char¬ 
acteristic into the model. Generally, n(k, rn) is a monotoni- 
cally increasing function of m for any given degree k, which 
characterizes the reinforcement effect through non-redundant 
information memory. If n(k, m) is a constant, no such re¬ 
inforcement effect exists. In this case, if the adopted state 
is regarded as the infected state in epidemiology, our model 


will reduce to the standard SIR (susceptible - infected - recov¬ 
ered) model fT 2 ;], where social reinforcement effect and non- 
Markovian properties are not present - a key difference be¬ 
tween biological and social contagions. Empirical researches 
indicate that the adopted individuals may lose interest in the 
behavior f34], which is also concerned in the binary social dy¬ 
namics [35], 3g]. At the same time step, we thus assume that 
each adopted individual loses interest in transmitting the be¬ 
havioral information and becomes recovered with probability 
7 . The spreading dynamics terminates once all adopted indi¬ 
viduals have become recovered. 

By setting the parameters, our stochastic model can gen¬ 
erate either Markovian or non-Markovian processes, thereby 
including a number of existing models on social contagions as 
different limiting cases. For example, if n(k, m) is a Heavi¬ 
side step function (i.e., if m is less than the adoption threshold 
T u , then 7 r is zero; otherwise, 7 r is unity), and setting A = 1.0 
and 7 = 0.0, we obtain the Watts threshold model [11]]. Once 
the thresholds of individuals and network topology are fixed, 
the cascade process in the Watts threshold model will be de¬ 
terministic, which is a trivial case of Markovian process. In 
addition, by choosing the dynamical parameters properly, we 
can map our model into some of the existing non-Markovian 
models. For instance, fixing A = 1.0 and letting n(k,m ) be 
a function of exactly one of the two quantities (i.e., adopted 
and susceptible individuals), we recover the synergy spread¬ 
ing model ll33ll . Similarly, if we allow n(k,m) to be a lin¬ 
ear fij or exponential f27] function of m and 7 = 1 . 0 , we 
can obtain distinct types of non-Markovian dynamics. Differ¬ 
ing from the models in Refs. |25i [271 in which each adopted 
individual only gets one chance to transmit the behavioral in¬ 
formation to every neighbor, in our model an adopted indi¬ 
vidual can try to transmit the information many times until 
he/she becomes recovered state or transmits the information 
successfully. 

In our study, we concentrate on the so-called spreading 
threshold model before turning to more generalized social 
contagion models. In the spreading threshold model, an indi¬ 
vidual u adopts the behavior only when the number of pieces 
of non-redundant information that u possessed exceeds the 
adoption threshold T u . This means that the adoption prob¬ 
ability n(k,m) is a Heaviside step function, which has the 
same form as in the Watts threshold model m. There are, 
however, key differences between the two types of threshold 
models. Firstly, differing from the Watts model in which the 
adoption threshold is the corresponding fraction of neighbor¬ 
ing nodes, the adoption threshold in our model is expressed 
in terms of the absolute number of neighboring nodes, as in 
bootstrap percolation |37|] and self-organized criticality mod¬ 
els HI. Secondly, in the Watts model, each individual can 
obtain information about the states of all its neighbors “in¬ 
stantaneously” at each time step, but in our model individuals 
are able to know the neighboring states only through trans¬ 
mission of the information. Thirdly, in the Watts model an 
individual is permanently interested in the behavior even af¬ 
ter its adoption, while we assume more realistically that in¬ 
dividuals having adopted certain behavior may lose interest 
in it and never spread the corresponding information, which 
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is quantified by the abandon probability 7. Note that if the 
threshold of Watts model is expressed as the absolute number 
of neighbors who have adopted the behavior, there will only 
exist the second and third differences. These three differences 
are consequences of introducing the non-redundant informa¬ 
tion memory characteristic into our model, better capturing 
the essential dynamics of social contagions in the real world. 


III. THEORY 

We first develop a unified edge-based compartmental the¬ 
ory to analyze our general social contagion model with re¬ 
inforcement mechanism based on non-redundant information 
memory characteristic. We then systematically investigate 
how the memory affects the social contagion process in a spe¬ 
cific model, the spreading threshold model. In this theory, 
we assume that the networks have large network sizes, sparse 
edges, and no degree-degree correlations, and the contagion 
dynamics evolves continuously. Mathematically, a contagion 
process can be described by three variables: S(t), A(t) and 
R(t), which are the densities of the susceptible, adopted, and 
recovered individuals at time t, respectively. The states of all 
individuals remain unchanged when t —> 00, and R( 00) is 
the final fraction of individuals that have adopted the social 
behavior. 


A. General theoretical framework 

Due to the non-redundant information memory characteris¬ 
tic, in a social contagion process there are strong dynamical 
correlations between the states of the adjacent nodes, making 
existing theoretical methods such as the mean-field theory @], 
percolation theory lITbh . and renewal process [ 30] inapplica¬ 
ble, especially for networks that are strongly structurally het¬ 
erogeneous. Using insights from Refs. II391 - I42I1 . we develop 
an edge-based compartmental theory to analyze social conta¬ 
gion dynamics in the presence of strong nodal state correla¬ 
tions. 

Let 9(t) be the probability that individual v has not trans¬ 
mitted the information to individual u along a randomly cho¬ 
sen edge by time t. In the spirit of the cavity theory 140:. 431, 
we disallow individual u to transmit any information to its 
neighbors but u can receive such information from its neigh¬ 
bors - u is in a cavity state. Initially, a fraction of po indi¬ 
viduals is in the adopted state, and none of them transmits 
the information to its neighbors, so 0 ( 0 ) = 1 for all edges. 
For simplicity in theory, we assume that the probability of not 
transmitting the information is identical for all edges, and dy¬ 
namical correlations doesn’t exist among neighbors of an in¬ 
dividual. At time t, a uniformly randomly chosen individual 
u of degree k in the cavity state has been exposed to m pieces 
of non-redundant information (i.e., u has received the infor¬ 
mation from distinct neighbors m times) with the probability 

tm(k, 9(t)) = (1 - po) ( * ) e(t) k ~ m [1 - (1) 


where the factor (1 — po) is the fraction of susceptible nodes 
initial. By time t, the susceptible individual u has received 
the information from m different neighbors. The probability 
that u has not adopted the behavior for time of receiving in¬ 
formation less than m is UJL 0 [1 — tt( k,j)\. Combining this 
factor and summing over all possible values of m, we obtain 
the probability that the individual u is still in the susceptible 
state at time t as 

k 

s(k,t)= Y 0m(M)n£L o [l - n(k, j)}. (2) 

m= 0 

Taking into account different degrees in the network, we ob¬ 
tain the fraction of susceptible individuals (i.e., the probability 
of a randomly chosen individual is in the susceptible state) at 
time t as 

OO 

S( i) = X>(*MM)- (3) 

k= 0 

Analogously, the fraction of individuals with m pieces of in¬ 
formation at time t is 

OO 

$(m,f) = y 'P(k)<f>m{k,0{t)). (4) 

k—0 

A neighbor of individual u may be in one of susceptible, 
adopted, or recovered states. We can thus further express 6(t) 
as 

8(t) = €s(t) + £ a (*) + £ r (£), ( 5 ) 

where £g(t) [£a(t) or £r(£)] is the probability that a neigh¬ 
bor of the individual u in the cavity state is in the susceptible 
(adopted or recovered) state and has not transmitted the infor¬ 
mation to individual u through an edge by time t. Note that 
the three quantities are unknown, which are to be solved. 

If a neighboring individual v of u is initially in the sus¬ 
ceptible state with probability 1 — po, it cannot transmit the 
information to u. Individual v can get the information from 
its other neighbors, since u is in a cavity state. At time t, the 
probability that individual v has received m pieces of non- 
redundant information is 

Tm(k',e(t)) = ~ i y{t) k '- m ~ l [\ - (6) 

where k 1 is the degree of v. Similar to Eq. ([J]), individual v 
will still stay in the susceptible state at time t with the proba¬ 
bility 

k'—l 

0(fc',0(f)) = Y T m (V,f)n™ 0 [l-7r(V,j)]. (7) 

m— 0 

For uncorrelated networks, the probability that one edge from 
individual u connects with an individual with degree k' is 
k'P{k)/{k), where (k) is the mean degree of the network. 
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Summing over all possible k', we obtain the probability that u 
connects to a susceptible individual by time t as 


&(*) 


(1 k ,k'P(k)Q(k',e(t)) 


( 8 ) 


The information spreading process as described in Sec.HIlsug- 
gests that two events need to occur in order for the growth of 
(1) with probability 1 — A an adopted neighbor has not 
transmitted the information to u via their connection and ( 2 ) 
with probability 7 the adopted neighbor has been recovered. 
Taking these into consideration, we get 


= 7 (1 - ( 9 ) 

At time £, the rate of change in the probability that a random 
edge has not transmitted the information is equal to the rate at 
which the adopted neighbors transmit the information to their 
susceptible neighboring individuals through edges. Thus, we 
get 


d,0(t) 

dt 


-A U(t). 


Combining Eqs. © and ( 1 1 Ob . we obtain 


( 10 ) 




7 [l- 6 >(f)](l~A) 

A 


( 11 ) 


Substituting Eqs. ® and (fill into Eq. ([5}, we get an expres¬ 
sion for in terms of 9{t). Doing so, we can rewrite 
Eq. (fTOl) as 


~dF - ~ X[m - (1 " Po) -<jfe)- 

+ 7[1 - - A). 


Note that the rate dA(t)/dt is equal to the rate at which S(t) 
decreases, because all the individuals moving out of the sus¬ 
ceptible state must move into the adopted state minus the rate 
at which adopted individuals become recovered. We have 


dA(t) 

dt 


dS(t) 

dt 


lA{t) 


( 13 ) 


and 


dR(t) 

dt 


7 Mt)- 


( 14 ) 


Equations O-© and (fiTTi-ifT-fli give us a complete and gen¬ 
eral description of social contagion dynamics, from which the 
density for each type of individual in each state at arbitrary 
time step can be calculated. 

Say we are interested in the steady state of the social con¬ 
tagion dynamics. Setting the right side of Eq. (ITU) to be zero, 
we get 


0(oo) = (1 - p 0 ) 


k'P{k')Q{k',9{ 00 )) 


(k) 


+ 


7[1 ~ 0(oo)] (1 - A) 


( 15 ) 


where 0(fc',0( 00 )) is a nonlinear function of 9( 00 ). Note 
that 9(t) decreases with t, as the individuals in the adopted 
state persistently transmit the information to their neighbors. 
Thus in simulations, only the maximum value of the stable 
fixed point (if there exist more than one stable fixed points) of 
Eq. (IT5l ) is physically meaningful. Substituting this value into 
Eqs. (Q])-([3]), we can obtain the value of the susceptible density 
£>( 00 ) and the final behavior adoption size R( 00 ). 

As in epidemic spreading, the condition under which out¬ 
break of behavior adoption occurs is of interest. Similar to 
analyzing epidemic spreading, we can obtain the critical con¬ 
dition by determining when a nontrivial solution of Eq. (1151) 
appears, which corresponds to the point at which the equation 

#(oo), POj T, 7,A] = (l-p 0 ) E *' fc,p ( fc g( fc '- e (°°)) 

+ tMMIBI -fl(oo) 

is tangent to horizontal axis at the critical value of 9 c (oo). 
The value of 9 C ( 00 ) denotes the critical probability that the 
information is not transmitted to u via an edge at the critical 
transmission probability when t —> 00 . This way we obtain 
the critical condition of the general social contagion model as 


dg 

d9{ 00 ) 


l(? c (oo) 0 * 


( 16 ) 


From Eq. (fl6l) . we can calculate the critical transmission prob¬ 
ability: 


A c 


7 

A + 7 — 1 


( 17 ) 


where 


A = (1 - p 0 ) 


V 1 T3tu\ dQ(fc',6>(oo)) I 
l^k' k ) dd(oo) 


(k) 


From Eqs. ©-([7]). we obtain the expression of as 


dQ{k\9{ 00 )) _ /fc' - 1\ 
d6{ 00 ) ^- 7 , \ to J 

x 7 m —0 x 7 

x {(*' - m - l)0(oo) fe, - m " 2 [l - 0(oo)] m 
- m9(oo) k '- m - 1 [l - f^oo)]™- 1 } 
xn“ 0 [i -n(k',j)\. 

(18) 

Numerically solving Eqs. (fl5l ) and (IT71) - (IT8l) . we can get the 
critical value of the transmission probability A c for any given 
adoption function w(k, m). We see that A c is correlated with 
the dynamical parameters such as the adoption probability 
7 t(/c,to), the initial seed size po an d the abandon probability 
7, as well as the topological parameters of the network [e.g., 
the degree distribution P(k ) and the mean degree (fc)]. 


B. Spreading threshold model 

We now apply the general theoretical framework developed 
in Sec. IIII Al to analyzing a specific class of social contagion 


A 
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model - spreading threshold model. In this model, the adop¬ 
tion function tt( k, to) is a Heaviside step function: 


7r(fc, to) 


1, TO > T k , 
0 , TO < T k , 


(19) 


where Tk is the adoption threshold of individuals of degree k. 
Here the adoption probability n(k, to) is only a function of k 
and m. Further investigations on general model, incorporat¬ 
ing individuals’ inherent characters such as age and habit, are 
called for. Utilizing Eq. ( flhl ). we can write Eqs. (O and Q as 

T k -1 

s(k,t)= ^2 <t>m{k,0(t)) ( 20 ) 

vn —0 


and 


Ty -1 

®(k',0(t)) = ]T r m {k',e(t)), (21) 

m= 0 

respectively. Similarly, Eq. ([13} becomes 

^ = A Um - (22) 

at 

where 

OO Tjfc — 1 / 7 \ 

V’W = (i-Po)53 p ( fc ) ( m ) 

k— 0 m—0 ^ ' (23) 

x {(k - m)e{t) k - m - 1 [i - e(t)] m 

-md{t) k ~ m [l -0(f)]™- 1 }. 

The critical condition can be determined using Eq. (fl 6 l) . For 
the simple case where the fraction of randomly chosen initial 
seeds is vanishingly small (i.e., po —> 0 ) and all individuals 
with different degrees have the same adoption threshold T, 
Eq. (H5l i has one trivial solution: 0(oo) = 1. At the critical 
point, the function g[6(oo), po, T, 7, A] is tangent to horizontal 
axis at 0(oo) = 1. For T = 1 , using Eqs. (fl5l)-(fT8l). we obtain 
the continuous critical transmission probability as 


11 = 7(fc) 

(fc 2 ) - 2 (k) + 7 <fc) ’ 


(24) 


which has the same form as the epidemic outbreak threshold. 
However, for T > 1, the function g[9(oo), po, T, 7 , A] can 
never be tangent to horizontal axis, suggesting that a vanish¬ 
ingly small fraction of initial seeds cannot cause a finite frac¬ 
tion of the individuals to adopt the behavior. 

Now suppose that po is not vanishingly small so that 
0(oo) = 1 is no longer a solution of Eq. ( IT5} . In this case, 
regardless of the values of other parameters, a finite frac¬ 
tion of individuals will adopt the behavior. It is thus rea¬ 
sonable to focus on how non-redundant information mem¬ 
ory characteristic affects the dependence of the final behav¬ 
ioral adoption size R( 00 ) on the transmission probability A, 
which can be obtained from the roots of Eq. ( fl5l ). We are 
particularly interested in finding out whether the dependence 




e 

FIG. 2. (Color online) Illustration of graphical solutions of 

Eq. (HB. For ER random networks, (a) continuously increasing be¬ 
havior of R(oo) with A for T = 1, (b) discontinuous change in 
R( 00 ) for T = 3. The black solid lines are the horizontal axis and 
the red dots denote the tangent points. Other parameters are po =0.1 
and 7 = 1.0. 


is continuous or discontinuous. Note that the number of 
roots Eq. (fl5l) is odd (including multiplicity) for any param¬ 
eters, because of the g[0( 00 ), po, T, 7 , A] <0 for 0(oo) = 1 
and g[0(oo), po, T, 7 , A] > 0 for 6 ( 00 ) = 0. As shown in 
Figs. [2a) and [2jb), numerical calculations indicate that the 
number of roots is either 1 or 3. When we fix all the parame¬ 
ters except A, if Eq. (115} has only one root for different values 
of A [Fig. Ha)], R( 00 ) will increase continuously with A. If 
the number of roots of Eq. ( 1 13 [ 1 depends on A, as shown in 
Fig.[2b), there will be three roots (fixed points), which means 
a saddle-node bifurcation occurs (44). The bifurcation analy¬ 
sis of Eq. C3 reveals that the system undergoes a cusp catas¬ 
trophe: Varying A, one finds that the physically meaningful 
stable solution of 6 ( 00 ) will suddenly jump to an alternate out¬ 
come. In this case, a discontinuous growth pattern of f?( 00 ) 
with A emerges, and the critical transmission probability X 1 ,, 
at which the discontinuity occurs can be obtained by solving 
Eqs. (HU) and (E}-®- 

The discontinuous behavior in R(po) versus A can be un¬ 
derstood by using a specific example, e.g., by setting T = 3. 
As shown in Fig. |2b), for different values of A, the func¬ 
tion g[0(oo), po, T, 7 , A] is tangent to horizontal axis at A^ ~ 
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0.5811. For A < A', if Eq. (TBT i has 3 fixed points then the 
solution will be given by the largest one (since only this value 
can be achieved physically). Otherwise, the solution is the 
only fixed point. For A = A 4 , the solution is given by the tan¬ 
gent point. For A > A 4 , the only fixed point is the solution of 
Eq. (IT5l ). In this case, the solution of Eq. < IT5l ) changes abruptly 
to a small value from a relatively large value at A = A 4 , lead¬ 
ing to a discontinuous change in R( oo). 

Finally, to determine the critical system parameter value of 
9 S ( oo), across which the dependence of R(oo) on A changes 
from being continuous (discontinuous) to discontinuous (con¬ 
tinuous), we can numerically solve Eqs. (IT5l) and ( 1 1 61 ) together 
with the condition 


d 2 g[9(oo), Po, T, 7, A] 
d0 2 (oo) 

From Eq. CHJ. we obtain 


= 0 . 


P o = 


F 


(25) 


(26) 


where F = J2k' k'P{k ') de • Using Eqs. (01 and 
<ED- we get 

dQ 2 (k',9{ oo)) _ Ya' (k' - 1 
d0 2 {oo) ~ V m 

x {( k' — TO — l)[(fc' — TO — 2)9(oo) k ~ m ~ 3 
x (1 — 6»(oo)) m - TO6»(oo) fe ' _m_2 (l - 9{oo)) r 

- m[{k' - to - l)0(oo) /c '- m " 2 (l - 9{oo)) m ~ 

- (to - l)0(oo) fe,_m_1 (l - 0(oo)) m - 2 ]}. 

(27) 

Combining Eqs. ( IT5ll . ( fl6l> and (l25l >. we get the value of 
9 S ( oo). For fixed T and P(k), the critical values of other 
system parameters e.g., A* and pg, can then be determined. 

IV. NUMERICAL VERIFICATION 
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I-IG. 3. (Color online) Spreading threshold model on ER net¬ 
works. (a) Average densities of susceptible, adopted, and recovered 
populations, denoted by S(t) (black squares), A(t) (blue up trian¬ 
gles), and R(t) (red circles), respectively, versus time, (b) Final be¬ 
havior adoption size R( oo) versus the transmission probability A for 
T = 1 (black squares), T = 2 (blue up triangles), T = 3 (red cir¬ 
cles), T = 4 (dark green diamonds), and T = 5 (light green stars) in 
the steady state, (c) Simulation results of NOI (number of iterations) 
as a function of A with T = 2 (blue solid line), T = 3 (red dashed 
line) and T = 4 (dark dash dotted green line). The error bars indicate 
the standard deviations. The lines in (a) and (b) are theoretical pre¬ 
dictions based on solutions of Eqs. 0-0 and 0D-CE3- In (a), we 
set A = 0.8, po = 0.1, T = 3, and 7 = 0.5 (so as to obtain longer 
evolution time), while in (b) and (c), we set po = 0.1 and 7 = 1.0. 


We perform extensive simulations on the spreading thresh¬ 
old model. In our simulations, we use Erdos-Renyi (ER) 
network model l|45i] and confi gura tion network model with 
power-law degree distribution [32j], where the network size 
and mean degree are N = 10 4 and ( k) = 10, respectively. At 
least 2 x 10 3 independent dynamical realizations on a fixed 
network are used to calculate the pertinent average values, 
which are further averaged over 100 network realizations. We 
separately discuss the effects of dynamical and topological pa¬ 
rameters. 


A. Effects of dynamical parameters 

To be illustrative, we use ER random networks ll45tl . We 
first calculate the time evolution of the population densities 


for A = 0.8, po = 0.1, T = 3, and 7 = 0.5, as shown in 
Fig. 0a), where we observe that the density of the suscep¬ 
tible (recovered) individuals decreases (increases) with time, 
and reaches some final value for large time. The density of 
the adopted individuals decreases initially (due to the fact that 
the number of individuals who newly adopt the behavior is 
less than that of individuals who become recovered), then in¬ 
creases and reaches a maximum value at f « 5. These results 
agree well with the predictions from our edge-based compart- 
mental theory [see lines in Fig. 0a)]. 

We next study the final behavior adoption size R( 00 ) as a 
function of the transmission probability A for different val¬ 
ues of the adoption threshold T at another value of 7 = 1.0. 
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FIG. 4. (Color online) The final fraction of individuals in the sub- 
critical state on ER networks. <f>(T — 1, oo) versus A for T — 2 
(black squares), T — 3 (blue up triangles), T = 4 (red circles) and 
T = 5 (dark green diamonds). The lines are theoretical predictions 
based on solutions of Eqs. 0. 0 and CCD-CEi- Other parameters 
are 7 = 1.0 and p 0 = 0.1. 


As shown in Fig. 0b), increasing T impedes individuals from 
adopting the behavior, since a larger value of T requires the 
individual to be exposed with more information from distinct 
neighbors to affirm the authority and legality of the behav¬ 
ior. As a result, individuals hardly adopt the behavior when 
the adoption threshold is relatively large (e.g., T > 5). Lines 
from the theory in Fig. 0b) are very consistent with these sim¬ 
ulation results. Through the bifurcation analysis of Eq. (ED), 
we note that the adoption threshold affects strongly the man¬ 
ner by which R(oo ) increases with A for T < 4. As shown 
in Fig. 0b), for some small adoption threshold (e.g., T = 1), 
R( 00 ) increases continuously with A. However, for a slightly 
larger adoption threshold (i.e., T > 1), the R( 00 ) versus A 
pattern becomes discontinuous, exhibiting an abrupt increase 
at some critical value X' r . The statistical errors are generally 
small except for A close to X', (for this reason and for figure 
clarity the error bars will not be shown for subsequent figures). 
The theoretical value of Xj. can be calculated from Eqs. (fl5l i 
and The critical value can also be estimated by ob¬ 

serving the number of iterations |4(| @3l (denoted as NOI, 
where only those iterations at which at least one individual 
adopts the behavior are taken into account). We observe that 
the NOI exhibits a maximum value at A^, e.g., for T = 2, 3, 
and 4, as shown in Fig. 0c). Overall, there is a remarkable 
agreement between theory and numerics in terms of the quan¬ 
tities R ( 00 ) and A'. Through extensive simulations and theo¬ 
retical predictions, we know the abandon probability doesn’t 
qualitatively affect the growth patterns of f?( 00 ), so it is set as 
7 = 1.0 in the rest of this paper. 

It is useful to identify the key factors that affect the depen¬ 
dence of R( 00 ) on A. To obtain an intuitive understanding 
of the phenomenon of abrupt increase in R( 00 ) as A passes 
through a critical point, we focus on the individuals in the 
subcritical state. An individual u in such a state has received 
the information but has not yet adopted the behavior, and the 



FIG. 5. (Color online) Dependence of final behavior adoption 
size on initial seed size and transmission probability. For spread¬ 
ing threshold model on ER random networks, color-coded values of 
R( 00 ) from numerical simulations (a) and theoretical solutions (b) 
in the parameter plane (po,A), where the theoretical values are from 
solutions of Eqs. 0-0 and The numerically obtained 

critical values of the transmission probability, A^ (white circles), are 
from the NOI method, and the corresponding theoretical values blue 
dashed line) are from Eqs. 0D> and 0-00. In each subfigure, 
three regions are shown: only a vanishingly small fraction of indi¬ 
viduals can be exposed to adopt the behavior in region I, in region 
II R(o o) grows discontinuously with A and a finite fraction of indi¬ 
viduals adopt the behavior above A^, and R( 00 ) grows continuously 
to a large value in region III. The vertical white solid lines and dash 
dotted yellow lines separate the plane into the three regions, which 
are predicted from our edge-based compartmental theory. Other pa¬ 
rameters are 7 = 1.0 and T = 3. 


number of information pieces from distinct neighbors is pre¬ 
cisely one less than the adoption threshold. Say at the time u 
has received information from his/her neighbors except neigh¬ 
bor v. Now assume that v has adopted the social behavior and 
transmits the information to u successfully so as to cause u to 
adopt the behavior. In turn, u will transmit the information to 
his/her susceptible neighbors with probability A, which will 
cause some subcritical state neighbors to adopt the behavior 
accordingly, and so on, potentially leading to an avalanche of 
behavior adoption. If the system has a relatively large num¬ 
ber of individuals in the subcritical state, a slight increase in 
the number of individuals who adopt the behavior, e.g., by 
increasing the value of A slightly, may cause a sudden and 
large number of such subcritical state individuals with infor- 
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mation pieces greater than their threshold, leading to a dis¬ 
continuous “jump” in the value of f?(oo). The above intuitive 
understanding is further proved by numerical simulations and 
theoretical predictions in Fig. [4] For 2 < T < 5, the final 
fraction of individuals in subcritical state <f>(T — l,oo) first 
increases with A below <f>(T — 1, oo) reaches a maximum 
at the A, 7 ; and a slight increment of A induces a finite fraction 
of $(T — 1, oo) to adopt the behavior simultaneously, which 
leads to a discontinuous jump in the value of R( oo). When 
this social reinforcement effect is not present [e.g., T = 1 in 
Fig- Ob)], there are essentially no individuals in the subcrit¬ 
ical state. In this case, R( oo) increases continuously with A. 
We mention that the mechanism underlying the discontinuous 
increase in f?( oo) in our spreading threshold model is similar 
to that responsible for phenomena such as explosive percola¬ 
tion [5U], bootstrap percolation j37tl. A-core percolation (49j 
and explosive synchronization [50]]. 

We further investigate the role of the initial seed size po 
in social contagion dynamics for relatively larger values of 
T, e.g., T = 3. As shown in Fig.0 we see that f?( oo) in¬ 
creases with p 0 , since individuals in the network have more 
chances to be exposed to the information. Based on the val¬ 
ues of R(oo), we can divide the phase diagram into local 
(po < 0.04) and global (0.04 < po < 1) behavior adoption 
regions, where in the former (i.e., region I), only a vanish¬ 
ingly small fraction of individuals can be exposed to adopt 
the behavior and, in the latter including regions II and III, a fi¬ 
nite fraction of individuals adopt the behavior and a crossover 
phenomenon occurs in the dependence of R( oo) on A. The 
crossover phenomenon means that the dependence of R( oo) 
on A can change from being discontinuous to being continu¬ 
ous. More specifically, the saddle-node bifurcation of Eq. ( TBT i 
occurs for 0.04 < p 0 < 0.15 (region II in Fig. [5}. thus R( oo) 
versus A is discontinuous; R( oo) versus A is continuous for 
0.15 < po < 1 (region III in Fig. [5}- as the saddle-node bifur¬ 
cation disappears. The crossover phenomenon originates from 
the fact that the number of individuals in the subcritical state 
decreases with po- At the crossover or switching point pg, 
as indicated by the vertical yellow dash dotted line in Fig. 0 
the behavior of R( oo) versus A changes from being discon¬ 
tinuous to continuous. The crossover point can be calculated 
analytically by solving Eqs. (IT5]i-(IT6t and (l25l i. We also find 
that A^ decreases with po, since a large value of po enhances 
the probability of individuals’ exposure to the information. In 
short, R( oo) versus the parameter plane (po, Aj_ shows a cusp 
catastrophe (i.e., the crossover phenomenon) [441 ■ Regardless 
of the size of the initial seeds, there is a good agreement be¬ 
tween numerically calculated and theoretically predicted be¬ 
haviors of R( oo). 


B. Effects of topological parameters 

We turn to elucidating the effect of network topological 
parameters on social contagion dynamics in our spreading 
threshold model. In fact, both the value of R( oo) and its pat¬ 
tern depend strongly on the mean degree and degree hetero¬ 
geneity of the network. To be concrete, we first examine ER 


random networks with different values of the mean degree (k), 
as shown in Fig.|6l where we see that R(oo) increases with ( k ) 
in general, since individuals with larger degrees have higher 
probabilities to be exposed to information from distinct neigh¬ 
bors, leading to a high likelihood that they adopt the behavior 
as well. By the bifurcation analysis of Eq. (IT5l ). we find that 
with the increase of (k), the growth pattern of R( oo) changes 
from being continuous to being discontinuous. For a small 
value of the mean degree (e.g., (k) = 5), only a small frac¬ 
tion of the individuals adopt the behavior, so R( oo) changes 
with A continuously. For a relatively larger value of the mean 
degree (e.g., (k) > 5), more individuals adopt the behavior, 
leading to a sudden, discontinuous increase in f?( oo) with A. 
As discussed in Sec. IIV A1 the “explosive” growth of R( oo) 
occurs whenever there is a finite but sizable fraction of indi¬ 
viduals in the subcritical state, which cannot happen when the 
mean degree of the network is small. We also observe that in¬ 
creasing the mean degree can reduce the value of the critical 
point A 7 , due to the corresponding increase in the number of 
individuals having relatively large degrees. 

We next study scale-free networks. Figure [7] shows, for 
T = 3, R( oo) versus A for (k) = 10. The uncorrelated 
networks are generated with the power-law degree distribu¬ 
tion P(k) ~ fc -1/ (z/ being the degree exponent) according to 
the procedure in Ref. 113211 . where the maximum degree is set 
as kmax ~ y/N. We find that increasing the heterogeneity of 
network structure (by using smaller values of the degree expo¬ 
nent) promotes (suppresses) R( oo) for small (large) values of 
A. This result can be qualitatively explained as follows ll42] l: 
From Eqs. Q-Q, we know that hubs adopt the behavior with 
large probability. With the increase of network heterogeneity 
(i.e., decreasing u), the network has a large number of individ¬ 
uals with very small degrees and more individuals with large 
degrees. For small values of A, more hubs for small v facil¬ 
itate behavior spreading as they are more likely to adopt the 
behavior. But for large values of A, a large number of indi¬ 
viduals with very small degrees have a small probability to 
adopt the behavior, resulting in smaller values of R.(oc). By 
the bifurcation analysis of Eq. ( fl5l) . we also observe that the 
system has a critical degree exponent v s « 4.0, below which 
R(oo) versus A is continuous but above which the variation 
is discontinuous. That is, as the network becomes more het¬ 
erogeneous, we expect a change in the dependence of R( oo) 
on A from being discontinuous to continuous, since the exis¬ 
tence of strong degree heterogeneity can not make individuals 
in the subcritical state adopt the behavior simultaneously. We 
also note that the critical point A) decreases as the network be¬ 
comes more heterogeneous. Again, there is a good agreement 
between theoretical and numerical results. 


V. ALTERNATIVE MODELS OF SOCIAL CONTAGION 
DYNAMICS 

The edge-based compartmental theory developed in SecIlHl 
can be applied to more general social contagion dynamics 
with reinforcement effect derived from non-redundant infor¬ 
mation memory characteristic. Here, we demonstrate the use 
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FIG. 6. (Color online) Effect of mean network degree on social 
contagion dynamics. For ER random networks, R( oo) versus the 
transmission probability A for mean degree ( k ) = 5 (gray squares), 
( k) — 8 (blue up triangles), ( k ) = 10 (red circles), ( k ) = 15 (dark 
green diamonds), and (k) = 20 (light green stars). Other parameters 
are po = 0.1, 7 = 1.0 and T = 3. The lines are theoretical values 
of R( 00 ) calculated from Eqs. |T]l-l[3} and d 1 2t )- d 1 4b . 
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FIG. 7. (Color online) Effect of network heterogeneity on social 
contagion dynamics. For scale-free networks, R( 00 ) versus A for 
degree exponent v = 2.2 (gray squares), v — 3.0 (blue up triangles), 
v = 4.0 (red circles) and v —> 00 (dark green diamonds). The 
case for v —> 00 reduces to a random regular network with identical 
degree. The lines are theoretical values of R( 00 ) calculated from 
Eqs. (Q}-i|3]i and GD-®. Other parameters are po = 0.1, 7 = 1.0 
and T = 3. 


of our theory in analyzing two alternative, somewhat more 
complicated social contagion models: ( 1 ) correlated spread¬ 
ing threshold model in which the adoption threshold of each 
individual is correlated with his/her degree and ( 2 ) a general¬ 
ized social contagion model in which the behavior adoption 
probability n(k, to) is a monotonically increasing function of 

TO. 


A. Correlated spreading threshold model 


In reality, whether an individual adopts certain social be¬ 
havior depends on his/her personal characters such as age and 
habit, which are reflected by the corresponding degree in the 
social network. As a result, there is typically some corre¬ 
lation between an individual’s degree and his/her ability to 
adopt new social behaviors triggered by crossing the adoption 
threshold. For simplicity, we use a recently introduced rela¬ 
tion 1 H to account for the correlation between individual i’s 
adoption threshold and degree hi, as 



Kmax 


(28) 


where k max is the maximum degree, A a and a are two ad¬ 
justable parameters. For a = 0, the adoption threshold is un¬ 
correlated with the degree, and all individuals in the network 
share the same adoption threshold. For a > 0, the adoption 
threshold is positively correlated with the degree, i.e., individ¬ 
uals with larger degrees have higher adoption thresholds, and 
the opposite occurs for a < 0 . 

To investigate the effects of varying a on social conta¬ 
gion dynamics using the spreading threshold model, we set 
the mean adoption threshold (somewhat arbitrarily) to be 
(T) = 3. The sum of the adoption threshold in the network is 
T s = Eti T i- For a = 0, we have T s = (T)N = A Q= 0 iV. 
Further, we get 


A a 


A a=0 Nk“, 

X^ N k a 
2~ii =1 N; 


(29) 


Evidence in terms of the quantity f?( 00 ), which supports our 
edge-based compartmental theory for varying threshold as 
given by Eq. (l28l >. is presented in Fig. [ 8 ] We observe a reason¬ 
able agreement between the theoretical predictions and simu¬ 
lation results. Note that a affects not only the value of R{ 00 ) 
but also its dependence on A. In particular, for a > 0, in¬ 
creasing a causes the critical point \' r first to increase then to 
decrease. This result can be qualitatively explained by not¬ 
ing that, slightly larger values of a (e.g, a = 1 ) can cause 
the individuals whose degrees are near the mean degree of the 
network to hold larger adoption threshold. However, much 
larger values of a (e.g., a — 2 ) will generate hubs with larger 
adoption threshold, thereby reducing the adoption threshold 
for the individuals with degrees near the mean degree. Since, 
in a random network, the degrees of most individuals are close 
to the mean degree, this causes the non-monotonic change in 
%■ 

For a < 0, decreasing a facilitates individuals’ adopting 
the behavior, and the dependence of f?( 00 ) on A changes from 
being discontinuous to continuous by the bifurcation analy¬ 
sis of Eq. (US. Decreasing a makes individuals with small 
(large) degrees to hold larger (smaller) adoption thresholds 
than the case of a = 0. As a result, the values of R{ 00 ) 
are smaller than those for a = 0 in the large A regime. Since 
individuals with small degrees have relatively large adoption 
threshold, they have more difficulty in adopting the social be¬ 
havior, further decreasing the number of individuals in the 
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FIG. 8 . (Color online) Effect of degree-correlated spreading 
threshold on social contagion dynamics. For ER random networks, 
final adoption size R(o o) versus the transmission probability A for 
a = —2 (gray squares), a = — 1 (blue up triangles), a = 0 (red 
circles), a = 1 (dark green diamonds), and a = 2 (yellow left trian¬ 
gles). Other parameters are (k) = 10, ( T ) = 3, and 7 = 1.0. The 
lines are theoretical values of R( 00 ) from solutions of Eqs. ([TJ (O 
and d 12b - J 14b with adoption threshold given by Eq. \29\ . 


subcritical state and making the discontinuous behavior in 
R(oo) to disappear. 


B. A generalized social contagion model 

Recently, Centola performed an interesting experiment of 
the health behavior spreading in an online social network, and 
found that the behavior adoption probability is a monotoni- 
cally increasing function of m Q], but not the trivial case of 
Heaviside step function in the spreading threshold model and 
Refs. H0]. Therefore, we assume that a susceptible individ¬ 
ual adopts the behavior with probability 

7t (fc, m) = 1 — (1 — e) m , (30) 

where m is the accumulated times that the individual has been 
exposed to different sources, i.e., he/she has received the in¬ 
formation m times from the distinct neighbors, and e is the 
unit adoption probability. We can also use the edge-based 
compartmental theory to analyze the dynamical process of this 
model by substituting Eq. (l30t into various equations that give 
the solutions of e.g., R( 00 ). In particular, we rewrite Eqs. (0 
and 0 as 

k 

s(k,t) = E 1 - (31) 

m =0 

and 

k'-l 

B(fc',0(f)) = E r m (k\6(t))(l - e )ESloJ, (32) 

m =0 




FIG. 9. (Color online) Results from a generalized social contagion 
model. For ER networks, R(o o) versus the unit adoption probability 
e for A = 0.3 (a) and the transmission probability A for e = 0.3 (b). 
Two values of po are used: po = 0.01 (gray squares) and po = 0.10 
(blue up triangles). Additional parameters are 7 = 1.0 and (k) = 10. 
In both panels, the lines represent the theoretical values of R( 00 ) 
obtained from solutions of Eqs. CIJ 0 J and G3-G3 w hh 4’{t) given 
by ([ir¬ 


respectively, whereas Eq. (IT3b has the same form as Eq. ( 1221 ). 
The different aspect is that we need to replace Eq. (l23l > with 

m = (i - po) E p ( fc ) E L) (! - e ) Er =° i 

fc=o m=o (33) 

x [(jfc - m)d(t) k - m - 1 [1 - 9(t)} m 

— md(t) k ~ m [l — 

Substituting Eqs. GD-G3 into the corresponding equations, 
we can obtain a theoretical understanding of the dynamical 
evolution of the generalized social contagion model. We ob¬ 
serve that R{ 00 ) varies with A continuously by the bifurca¬ 
tion analysis of Eq. (fl5l ). The theoretical values of R( 00 ) so 
predicted agree well with the simulated results, as shown in 

Fig. [9] 

VL CONCLUSIONS 

In social contagion dynamics, memory of non-redundant 
information can have a significant impact on the reinforce- 
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ment mechanism required for behavioral adoption. In par¬ 
ticular, the non-redundant information memory has two fea¬ 
tures: (1) repetitive information transmission on every edge is 
forbidden, (2) every individual can remember the cumulative 
pieces of non-redundant information. Social reinforcement in¬ 
corporating the memory characteristic is essential to describ¬ 
ing and understanding social contagions in the real world. In 
this paper, we first proposed a general social contagion model 
with reinforcement derived from this memory characteristic. 
Mathematically, a model based on such characteristic is nec¬ 
essarily non-Markovian. Previous works pointed out the diffi¬ 
culty to develop an accurate theoretical framework to analyze 
social contagion dynamics with only memory effect @], let 
alone models with non-redundant information memory char¬ 
acteristic. To meet this challenge, in this paper we developed 
a unified edge-based compartmental theory to analyze social 
contagion dynamics with non-redundant information memory 
characteristic. The validity of our theory is established by 
testing it using different social contagion models of varying 
complexity, different model networks. 

Through a detailed study of a comparatively simple model, 
the spreading threshold model, the effects of non-redundant 
information memory characteristic on the social contagion dy¬ 
namics can be quantified by the final adoption size R( oo) and 
its dependence on key parameters such as A. Especially, de¬ 
creasing the adoption threshold, increasing the initial seed size 
or increasing the mean degree of the network can facilitate 
adoption of social behaviors at the individual level, making 
the system less resilient to social contagions. The effect of 
structural heterogeneity on R( oo) turns out to be more com¬ 
plex in that, while making the network more heterogeneous 
can promote the spreading process, it impedes spreading for 
relatively large values of A. A striking phenomenon is that 
f?(oo) as a function of A can exhibit two characteristically dif¬ 
ferent types of patterns: continuous variation or sudden, dis¬ 
continuous changes, and a transition between the two patterns 
can be induced by adjusting parameters such as individuals’ 
adoption threshold, the initial seed size or the structural het¬ 
erogeneity of the network. For example, in order to change the 


dependence of R( oo) on A from being discontinuous to con¬ 
tinuous, we can decrease the individuals’ adoption threshold, 
increase the initial seed size or make the network more het¬ 
erogeneous. We also find that the discontinuous pattern disap¬ 
pears when there is negative correlation between individual’s 
adoption threshold and his/her degree. The above crossover 
phenomena can be understood through the bifurcation anal¬ 
ysis in theory, and also justified by analyzing the subcritical 
individuals in simulations. 

To study social contagion dynamics in human populations 
is an extremely challenging problem with broad implications 
and interest. Our main contribution is a treatment of the non- 
redundant information memory characteristic that is intrinsic 
to real world dynamics of social contagions. Our unified edge- 
based compartmental theory gives reasonable understanding 
of the roles of the memory characteristic in shaping the 
spreading dynamics, which can be applied to analyzing dif¬ 
ferent dynamical processes such as information diffusion on 
computer networks. However, many challenges remain, such 
as incorporation of correlations between local structures (e.g., 
communities and motifs) into social reinforcement effect at 
the individual level, the impacts of redundant versus non- 
redundant information transmission, and further development 
of analytic methods to treat non-Markovian social contagion 
model on more realistic networks such as clustered 152. 531. 
multiplex ll54l - l58ll . and temporal networks ||59i[60l l ). 
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